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^ ■ Abstract 

We review and extend a technique for recovering a smooth function from its averages over 
P^P ' a wide class of curves in a general region of Euclidean space. The method is based on 

complexification of the underlying vector fields defining the transport and recasting the 
problem in terms of complex-analytic function theory. Conditions on the validity of prior 
formulae appearing in [HB10] as well as stability estimates are then discussed first for the 
■ case of vector fields with polynomial coefficients and later for more general cases. 

Keywords: X-ray transforms, explicit inversion, complex analysis, transport equation, har- 
monic calculus 



1 Background and Motivation 

The following filtered backprojection formula appeared in |Bal05] : 
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where f(z) is a smooth enough function and If its geodesic ray transform in the Poincare 
disc. Formula (PQ), once obtained, subsequently gave a holomorphic integrating factor to 
derive a similar, though more complicated, formula for the attenuated radon transform 
(AtRT) I a f on the same space. The method used in that paper, which we call the method of 
complexification, was an extension of one first used in [Nov02j and rests on the introduction of 
a complex parameter A into the governing transport equation and a subsequent analysis of the 
behavior of solutions in terms of this new parameter. 

Recently, in (HBlOj . we obtained a strikingly similar result to (JT|), namely 

f & = ^ £ r P(K0)X^H(I e f)(s(ze- ie ),e* e )d9 (2) 

The aim of this paper is to outline, briefly, the method which resulted in the above formula 
as well as to further classify the breadth of its validity. Although the original impetus for the 
above was the inversion of the AtRT, formula t^j is interesting in its own right. Much of the 
material in this article may be found in more detail in |Hoell] 

The structure of this paper is as follows. Section [2] contains the entire cast of characters 
needed throughout the article. Since our goal is partly expository, section [3] is a concise review 
of the method of complexification as well as a discussion of what we term "H-ness" and its 
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limitations. In section [5] we examine H-ness in more detail first for the case of polynomials 
and eventually for more general real-analytic vector fields. We generalize the situation in 
the penultimate section [7] where we present results relating H-ness to reasonable frequency 
constraints. 



2 Generic Preliminaries 

Our setup will be as follows. Let 7 : M. 2 3 (t, s) 1 — >- *y(t, s) G E C M 2 be a real- analytic 
diffeomorphism on a simply-connected domain E generating the linear, stationary transport 
operator 

d d 

X x = a(x,y) — + b(x,y) — , x=(x,y)G£ (3) 

We regard R 2 = C via the standard isomorphism so that 7 is identified with 7 1 (t, s) +i , y 2 (t, s). 
Defining complex w = j(t,s), we see that (w,w) are now (independent) complex coordinates 
on £. The regularity of the curves 7(t, s) show us that 7*^ is a non-degenerate field on E, 
X\ w = fi(w)-J^- + ^( w }jL where (0*X)(/) = X(<fi*f) defines the pushforward 0*. 

The equation of interest is the stationary transport boundary value problem X\ w u(w) = 
f(w), for w G £, f(w) G C%°(E) with lim^-^ u(w(t, s)) = 0, i.e. the BVP 

M^)^ + ^H^- = /H ; »eE (4) 

u \ d _ s = (5) 

One key difference between this formulation of the problem and that considered in |Bal05] is 
that there is, a priori, no immediately obvious object to "complexify" since we no longer have 
a parameterization of the points of impact on <9£ as was the case then. To circumvent the 
aforementioned difficulty, we appeal to the Riemann mapping theorem QNeh52[ IGK06J) which 
guarantees a unique biholomorphism z : E — > D + satisfying z(Q = 0, z'(Q > for £ G D + , 
where D + is the unit disc {z G C; \z\ < 1}. Since the Riemann mapping is conformal it is 
necessarily (infinitesimally) factorable (as in e.g. |Pal04] ) with respect to the subvarieties of 
integral curves of X z . Because of this equivalence between our initial domain E and the unit 
disc all further results will henceforth be presented in the disc. 

Since ~y*z maps M 2 into D + we use (z, z) as our coordinates on D + and have a new vector 
field X\ z = z*X\ z ( w ) where \x i-> {z*/j} o z~ l and likewise for \x. By a forgivable abuse of 
notation we denote {z*yu}^ o z^ 1 by fi(z) and {z*j2}j^r o z^ 1 by p,(z) so that 

x \z = Kz)-^ + fi(z)-^, zeD + , \fi\>o 

is our governing differential operator. Defining t(z) = z*w*t and s(z) = z^w^s, smooth func- 
tions on D + , the method of characteristics gives the solution to the BVP X\ z u(z) = f(z), 
u(z(— 00, s)) = as 

u{z) = {D 1 f){z) = ~ / f{z{t , s))sign{t{z) - t )dt (6) 

and since 7*z* = (z o 7)*. we define the X-ray transform of a function /(z) over the integral 
curves of X\ z , indexed by the transverse parameter s, to be 

(!/)(*)= / f(z(t,s))dt (7) 
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The main players we need at our disposal £1X6 8bS follows; 
Symmetric Beam Transform Q 

(D e i/>){z) = \ I ^(e id z(t , s(ze- i6 )))sign(t(ze' ie ) - t )dt , V G L 1 ( J D+) 
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Ray Transform 

m(s,e M ) = (I e iP)(s) = [ ^ e z(t,s))dt, 4> G L\D + ) 

JR 

Hilbert Transform 

(Hijj)(x) = -p.v. f ^-dy, V e L p (R), p > 1 

JrX-V 

as well as the classical Poisson kernel P(z, 9) = ijzpli^ra which for z G -D + , 9 E T generates the 
harmonic functions on the unit disc. Occasionally we will use the nonstandard notation Iq f 
to indicate the line integral of / over the integral curves of the vector field X. 

3 Complexification in a Nutshell 

The main result of [HB10] is the following. 

Theorem 3.1 Under suitable conditions on X\ and s(z) there exists a function Xi(z) on D + 
such that 

f(z) = ^ J * P(\,0)XtH(! 9 f){s{ze-*),e»)dB (8) 

provides a reconstruction for the function f based on the data Ief of ray transforms of f over 
the integral curves of Xq = 9*(X\ Z ). 

We now review the method outlined in that paper which was used to obtain this result. 
3.1 Symmetrizing and Symmetry-Breaking 

Let A = 9 E T = dD + and define the conformal map A : (z,z) — > (Xz,\z). If, for each s, 
$(•, s) is a set of integral curves of D + , then z -1 (A*$(-, s)) are conformally related curves in S. 
For A G D + /{0} we then consider X x = \*X\ Z to be the so-called "complexification of X\ z ". 
Explicitly, \*X\ Z takes the form Xz)X-^ + Xz)j-§= or X x = £(z, A)J^ + p(z, A) J| with 
j^(z, A) = n(z, A) = X*fi(z) and Xp(z, A) = fl(z, A) = X*p(z). 

Define X^ = ±i(—£(z, A) J^- + p(z, A) J|) as a vector field orthogonal to X\ when A = e %e . 
Namely, X e ■ Xq = ±i(\\(z, e i0 )\ 2 - \p(z, e i0 )\ 2 ) = in the standard inner product • : C 2 -> C. 
The prefactor of i makes Xjj~u(z, 9) real- valued and the choice of ± is determined by whichever 
satisfies the condition X^s > 0. Because X± = a(z)z*-^ for some real- valued this 
determines X^ uniquely and since we could just as well reparameterize with — s we will, without 
losing generality, avoid keeping track of signs by assuming X^ = z(— A)J^ + p(z, A)J|). 
Likewise, define s(z, A) and t(z, A) respectively as X*s(z) and X*t(z) for A G D + /{0}. □ 



1 N.B. e 10 and 9 will be used interchangeably, their meaning clear from context 

2 For functions k(z, A), ^| and k z are equivalent, as are ^| and k s , and we will use them interchangably. 
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3.2 Analysis, Asymptotics, and H-ness 

Our complexified transport equation now reads as follows; 

X x u(z,X) = f(z), A G D + (9) 

where it should be stressed that the parameter A has no obvious relation to the original particle 
transport that started this rigmarole. The method used to obtain Theorem 13.11 involves solving 
equation (jUj) for u(z, A) and showing analytic dependence of the solution on this parameter, i.e. 
d x u(z, A) = 0. 

A restricted class of vector fields known as type H was identified which ensure the ensuing 
steps work out as we need. The following is a revised version of that definition better suited to 
the purposes of this paper. H 

Definition A real vector field X\ z , complexified in the manner above 

KX\ Z = a(z, A)— + b(z, A) — , AG D+/{0} 
az az 

is said to be of type H if the following holds: 

1. a(z, A) is a holomorphic function of A for A G D + and has at least one zero A = \{z) G D + 

2. b(z, A) is a meromorphic function of A for A G D + and has no zeroes in D + /{0} 

3. f|j^y is a holomorphic function of A G D + 

4. s(z, A), d^' A ) ; meromorphic functions of A for A G D + 

where, as in the above, s(z, A) = \*s(z) is the complexified parameter specifying transverse 
foliation of the integral curves of X\. 

This condition H is the "suitable condition" mentioned in 13.11 and we will assume our vector 
field is of this type (i.e. a = [i and b = p). The \{z) mentioned previously are the zeroes of 
the complexified coefficient of our initial field. Note that condition H is strong insofar as 
holomorphy itself is a rather stringent condition. The above criteria will heretofore be called 
"H-ness". ' 

3.3 A Proof Sketch 

We give a scandalously brief sketch of the proof leading to (|SJ), highlighting where H-ness comes 
into play. First of all, by the third condition in 13.21 of nondegeneracy we see that the Jacobian 
ds(z) = \s z (z, A)| 2 — \sz(z, A)| 2 7^ holds on A G D + /{0} since the inequality 

guarantees that \s z (z, A)| 2 ^ on that same region. 

3 Notice that the maximum principle causes the third condition to follow automatically if, in the second 
condition, D + /{0} is replaced by D + which was the situation examined in |HB10| . Since the first and last of 
the above are in agreement with those already considered, we need only check that the middle conditions square 
with what we need. 
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We may therefore make a change of variables in s to get s*X\ = s*X\s(z, A)J|, whereby 
our fundamental equation X\G\(z; zq) = S(z — z ) is solved explicitly by 



d(t,s) 
d(z,z) 

G\(z, Zq) = —X 



, . (10) 
ir(s(z) - s(zo)) 

Checking against a bump function extends this to hold weakly at A — > and density shows 
that u(z, A) is holomorphic in A as needed. A similar argument works on u z (z, A) and u s {z, A) 
by invoking the final condition of H-ness. From here, Hilbert's relations on the boundary values 
of complex-analytic functions become viable as the following result shows. 

Proposition 3.2 

u ± (z,e* e )= lim u(z,\) = T^(HIef)(s(e- ie z),9) + (D e f)(z) 
where the Hilbert transform H is taken with respect to the first variable. 

The proof of this comes from an explicit analysis of [TO] with A = 1 — e and deriving the relation 









A=l 



AW(z, 1) = - JL^L X f 8 (z, 1) (11) 

It can be shown that by invoking the third condition of 13.21 we have 
sign(is'(z, 1) - is'(z , 1)) = sign(t(z, 1) - t(z Q , 1)) 
from which the Sokhotskyi-Plemelj formula allows us to obtain the advertised proposition. 

3.4 Reconstruction 

We use Proposition [321 together with the classical representation of complex- analytic functions 
on the unit disc. By definition of A« and Xjj~, one has iX\.u(z,\i) = Xyu(z, Aj) so that on 
equating real and imaginary parts we have 



^ P{K0)Xi{D e f){z)d9=-±- ^ P(X i ,9)X e H(Ief)(s(ze- id ),e i6 )d9 

27T J 47T J Q 



and 

1 r 2lT 



f(z) = —J P(X i ,9)X^H(I e f)(s(ze- M ),e^)d9 (12) 

which is the result we sought. Notice that H always denotes the Hilbert transform with respect 
to the s variable. 



4 Statement of Results 

The goal of this paper is to establish, over the next two sections, the following result (viz. 
Theorem 17.41 and Corollary 16.41 respectively) . 
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Theorem 4.1 Let e > be given. Suppose that fi(z,z) = Yl P +q>o a pq zP ^ 9 ^ s a real- analytic 
function on D + and that f G C%°(D + ). Define c r (z,z) = Y^ q - P = r a P q zP z q ■ Let l(z) and k(z) be 
the max and min respectively of the j such that Cj(z) ^ 0. Suppose that \*s, X*s z , and A*s g are 
meromorphic for A G D + . 

• If there are only finitely many nonzero Cj(z) , and both l(z) + k(z) + 2 > andO < |c&(z)| < 
|c/(z)| holds for all nonzero z G D + , then there exists a vector field Yq such that we have 
a perfect reconstruction 

»2tt 



M = ^ JJ P(Xt,e)Y^H(Ief)(s(ze- i6 ),e ie )d9 



If, for all z G D + /{0}, there exist infinitely many j G Z such that Cj(z) and Cj+i(z) are 

J j-S>00 I Cj (2) 



both nonzero, and if limsup,-_ Kjo | Cj+ ffi | < 1 then there exists a vector field Yq such that 



\\I Y 'f-I X f\\ LHD+) <Ce l<g<oc 

If, in addition, we have the Frechet bound \ \HsI x f — HIgf\\s < S(e) then there exists a 
function Aj : z — > \i(z) satisfying 



sup 

2S-D+ 



' < C(S(e)) (13) 



/(*) ~-^J ^ P(Xhe)Yi ± H(I e f)(s(ze- M ) 1 e i9 )de 



where HIgf is the Hilbert transform in the s variable of the trace of f over the integral 
curves of (/id + p,d) 

Clearly (120]) generalizes jSJ) in the sense that it allows an approximate reconstruction, to 
arbitrary accuracy, for a large class of vector fields. We will, in the sequel, prove Theorem 14.11 
and help explain just how broad its applicability is. The case of polynomial fields is as good as 
one could hope for. Our methodology is to establish results first for the case of polynomial vector 
field coefficients and later to reinterpret the terms Cj(z) as the frequencies of the complexified 
vector field's coefficients in the case of non-polynomial fields. 



5 Polynomial Vector Fields 
5.1 The polynomial space T(Q) 

Consider the case in which /jl(z) is a nonvanishing polynomial, i.e. fi(z, z) = J2 P + q >o a P q zP z q for 
z G f2 D {0}. The complexified coefficients of the field X\ around A = are then 

t{z,\) = J2 b pq (z)X^+\ P (z,X)= dp q (z)X p ~ q ~ 1 

p+q>0 P+Q>0 

with b pq = a pq z p z q and d pq = b pq . In order for H-ness to holcQ we will need that d\^{z, X) = 0, 
which a priori we do not have since q — p + 1 may very well be negative. If g — p + 1 > for all 
(p, g)-pairs then we are (provided we have roots and the rest of condition H) in the position 

4 For the moment we will be ignoring any possible problems with s(z, A) 
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of the previous section. If not, i.e. if p > q + 1 holds for some (p, q)-pair, then we proceed as 
follows. 

First of all, we will mostly be using the local irreducible exponents k, I given by 



k(z) = min (j) and l(z) = max (j) (14) 



where 



c r (z,z) = 



P,1 
q—p=r 



so that the Laurent expansion of fx(z, A) around A = is given by Y^! r =k c r{ z -> '■ Occasionaly 
we will need the global exponents defined as 

kn = = min k(z) and L = £(//) = maxZ(z) (15) 



Obviously —A" < and similarly < l(fj) < Nq, and our previous assumption is 

equivalent to the condition k(fi) + 1 < 0. Notice k > k^ and that I < depending on z G Q. 
To be clear, if there is no /i we are referring to the local irreducible exponents. Since \fi\ > we 
can be certain that k, I always exist (even if they may be equal). Also, k(0) = 1(0) = with 
co(0) = a 00 . 

We define the following polynomial space; 

r(fi) = {fi = ^2a pg z p z q ; \n\ > and \k(z)\ < I + 2} 
p+q 

for reasons which will be made clearer in the sequel. 
5.2 The rescaling scalar 

Consider the function w(z) = 2 — z^^ 1 — z^ -1 , which has two important properties: 

1. w(z) elforzeC 

2. < \w\ < 2 << oo for z G D + 

The first of the above guarantees that the field Y\ z = X\ z = a(z)^ + b(z)^= has the 
same integral curves as X\ . The second fact ensures that this rescaling introduces no artificial 
degeneracies into the field, in the sense that \a\ = |^| > 0. This amounts to a change in 
variables generated via 

dp 



t(t,s) 



o 2 — z\ k <*\ 1 (p,s) — z\ k v\ 1 (p,s) 



5 If < then we should use the complementary complexification A : (z,z) — > (jz,\z), for A € D + and 
get a holomorphic jA*/i. Since this is a situation which was dealt with in the previous section we may assume 
that l(n) > 0. 
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5.3 The first three conditions of H-ness 

Our first result towards establishing H-ness in the case of vector fields with polynomial coeffi- 
cients is the following simple lemma. 

Lemma 5.1 a(z, A) = X*a(z) is holomorphic for A G D + 

Proof We have a(z, A) = 2X \ k ,, \ -i J^^t- i_gdM-i) A 2(|fc M i-i) w here the numerator contains only 

\ k (0) + 1 

positive powers of X. Notice a(0, A) = — 2 ~ A. The quadratic formula shows that away 

from the origin z = the solution to w(z,A) = is given by |A(z)| = | — - — — r 

By the triangle inequality with z ^ we see that 



^IM- 1 



1<|I|2_|I^1_U|2|.< 

z z z 

with equality holding only when \z\ = 1, and therefore A ^ D + for \z\ < 1, and ipso facto 
w(z,X) 7^ for A G D + . 

We can now obtain a positive answer on the first criterion of H-ness. 
Proposition 5.2 // fi(z, z) G T(D + ) and if 

»2tt I 



1± ^/l 
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2 |2 ^ | V I "I \2 



log| C& (^)| <±-[ \og\Y,Cr{z)e i9 ^\ d9 
2n Jo ~t 



j=k 

for z G D + /{0} then the coefficient a(z,X) = 2A , fcp ] - 1 jgF^S^piJ^ i - d has a root M z ) G 
D + and t/ie /irst condition of H-ness is met. 

Proof Recall Jensen's formula for a meromorphic function h(z) with roots a v and (3 T in a 
region R = {z, \z\ < R}, 



log I MO) I = i f * log |Zi(Jte«)| d6 + V log %1 - V lo 6 



R 



(16) 



provided |/i(0)| 7^ 0, 00 fsee e.g. W. 64 , Nev7(tf ). Since the polynomial P/_jt(A) = c^^A^^ + 
Cfc + i(;z)A fc ~ fcfl+1 + • • • + Q(,2)A'~ fcM /ias no poles and since X = is not a root when k = k^, we 
may apply Jensen's formula in that case to P/_k(A) and R = 1 to yield 

1 /*2"7T j 

log|c fe (^)| = — / log| Vc^K^'-^l d^ + VloglA^)! 
27r J o j=k i 

where Pi~k{Xi(z)) = and the result is immediate. If k > k^, X = is a root of local order 
kn(z) — k^ and there's nothing to prove. At z = there is likewise nothing to prove. 

The next theorem uses similar arguments to address the second condition of H-ness. 
Proposition 5.3 If fi(z, z) G T(D + ) and if 

f\g|Vc^)e^M-;-2)| d9 < \og\- Cl {z)\ 
271 Jo 3=1 

then b(z, A) = 2 a 7m - 1^ 1 i 1 j I m - rj^p^rj " nonvamshmg for (z, A) G D + x D + / {0}. 



Proof Since fx was given as a polynomial we are guaranteed meromorphy of the term b(z, A). 
Looking at b(z, A) = ^^rz^^E^^T^m^) we see that b(z, A = 0) is nonzero for z ^ 
and \kfil —1 = 2 since q ^ 0. If z ^ and \k^\ — I < 2 then of course limm\o b(z, A) = oo 
with local order \k^\ — l^z) — 2. If z = then b(z, A) ~ v near |A| = 0. T/ie /acA; o/ vanishing 
of the denominator together with the way fi was complexified ensure that b(z, A) has no other 
singularities within D + x D + . 

In the Jensen formula (Tty) . if h(z) had a zero of order m at z = then h = h§z m + • ■ ■ in 
a vicinity of the origin. In that case the function T(z) = R has the same modulus on dR 
but is nonvanishing at the origin, its value there being R m ho. The Jensen formula applied to 
T(z) would yield 

log \ho\ = ^ log \h(Re i6 ) | dB + 5>g 1^1 - l0 S ~ m log R 
Assuming that there exists at least one Aj(z) G D + /{0} such that b(z,Aj) = 0, we use (17) 



with h(X) = ci(z)X^- l ~ 2 + ■■■ + c k {z)\\ k »\- k - 2 to get 

log|c,(z)| = — / 27r iog|Q(^)e*l-'- 2 ) + --- + c fe (z)e*l- fc - 2 )| d# + Vlog|A 
2n Jo a 



■3\ Z ) 



Whence 

log|c,(z)| < ±- [ "log | y^c 3 {z)e w ^-^\ d9 (17) 
T/ie result follows from the above inequality by contradiction. 

Remark Since f|j^y = ^jr r= _ fc ^*^_ 2 ? the origin A = is the only spot where analyticity may 

fail. But since G r(fl) we see that I + k + 2 > and |^y| < 1 keep | bounded as A — > so 
that by Riemann's theorem | is analytic on all D + as required in condition H. 

The preceding results combine in the following important corollary. 

Corollary 5.4 If fi(z,z) G T(D + ) and if for z ^ we have 

l0g|Cfc(*)| < log I Cj(2) I 

then X\ z = a(z, A)J^ + b(z, A) J| meets the first three conditions of H-ness. 

Since polynomials are the building blocks of real-analytic functions, we extend these results 
in the subsequent sections. 

6 H-ness in the space TLk,i{ty 

We begin with a definition extending our previous notion of r(fl). 

Definition A real-analytic function fi(z) = J2p+i=j a pq zP z q on Q D {0} with aoo 7^ belongs 

i>o 

to the space {?fc M ,z M (^) when 
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2. Both k(n) and denned as in ffTBT) are finite with < —1 and Z(/i) > 

3. — Ar(^f) < l(z) + 2 for all z G 

4. < |cfc(z)| < |q(z)| holds for all z G ^/{O} where c r (z, z) = J2 q - P =r a pq zP z 9 with k, I the 
local irreducible exponents of f }T4l 

We drop subscripts on Qk^i^iS^) since the notation Qk,i(P) is more concise and the global 
meaning is obvious. Clearly |cfc(0)| = |q(0)| = |ooo| > 0. The condition on absolute summabil- 
ity ensures that c r (z) is well-defined. The conditions guarantee we are left with a complexified 
\fi(zj, zX) which has a finite Laurent series in A. We have thereby established our main result 
with regard to polynomials. 

Theorem 6.1 Let 11 G Gk,i(D + ) and let w(z) = 2- z^^" 1 - z^"^ 1 . Then the vector field 

X X = X (V^LL— + ^) 9 \ 

w(z) dz w(z) dz 

satisfies the first three conditions of H-ness. 

With that in mind, we make the following 
Definition Denoting the meromorphic functions in A G £1 as we define 

H k ,i(n) = {fi G g kil (Q); s(z, A), ds(z, A), ds(z, A) G M(£l)} 

6.1 The Fourth Condition 

We now address the fourth and final condition of H-ness, namely meromorphy of s(z, A) and 
its z and z derivatives. To start with, this condition is already more relaxed than the initial 
three since meromorphy itself is less restrictive than holomorphy and there is no constraint on 
existence (or lack thereof) of roots. Secondly, for the space HL p (G, Q) defined as all / satisfying 
both 

1- g L p (n), Vz e G 

2. \\f\\ H L p{ G,u) = ll/IU^o) + suplliM^IU^);^ g G} < 00 
we have the result ( [Ren89j . Thm. 3.2) 

Theorem 6.2 For domain Q, Zq G Q and 11 G HL p (G,Q) for p > 2, \fi\ < 1 and V 'z G Q we 
have that if u(z) solves du = —fi(z)du on Q and u(z) has a zero/pole of order m at the point 
zq then 

u{z) = c{(z- zq) + b(z-^z- )} m + ck{(z - zo) + b(^~o)} m + 0(\z- z \ ±m+a ) 
for some a > and the ± picked according to whether zq is a zero or pole respectively. 
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The point is that at least locally we see an expansion for which, with A G D + constant (and hence 
||| 7^ 1), one should have meromorphy of the solution s(z, A) to A*(^yjj + ^) = 0- 

The above would constitute an expansion of r(z) = s(z,X) — s(zq,X) since that clearly has a 
zero at zq, although the order is not known a priori. Thus, for many non-pathological cases (i.e. 
excluding essential singularities, etc) r(z) would have the proper local expansion at all points 
in Q to satisfy meromorphy in A. Meromorphy of the derivatives would then follow. While this 
does not constitute a proof that condition 4 of H-ness is necessarily satisfied, it does constitute 
a proof modulo pathological cases. Clearly then Gk,i{&) C 'Hk,i{fy an d we can be sure that 
Hk,i(ty/Gk,i(ty is not too large. We now state the main result of this section. 

Theorem 6.3 If fj, — J2a pq z p z q G Hk,i(D + ), then the field 

X X = X (tQ— + & 9 ) 

w(z) dz w(z) dz 

with w(z) = 2 — z^^ 1 — z^vh 1 satisfies condition H 

The above result allows us to reconstruct functions over what are initially non-type H fields 
as in the following easy corollary. 

Corollary 6.4 If ^ e H k>l {D+), X^\^ z) = and = i9,(-$$& + and f G 

C °°(D + ) then 

f(z) = ^ jf" P(X l ,9)X^H(I e f)(s(ze- ie ),e w )d9 

where Igf is the ray transform of f over the integral curves ofYg = 6*(fid + fid). 

Proof Consider the equation X\ z u{z) = g{z) for g = G Cq°(D + ). Then by Ii6.3\) A* X\ z 
is type H and has zeros \{z) G D + . Thus, by (TJj) 

9{z) = ^ P(X i ,9)X^H(I e g)(s(ze- w ),e w )d9 

where Igg is the trace of g over the integral curves of Xg. However f was arbitrary in C^°(D + ) 
and since efw G M is both finite and nonvanishing on D + the integral curves of e~f(fid + p,d) 
and of X e are the same. In particular, under a change of variables, Igg = Igf. The result 
follows since s was unchanged. 

7 Some Harmonic Analysis: Onward and Upward 
7.1 The Projection Operator 

The Fourier expansion of a smooth function a(z, e 10 ) on the unit disc given is by 

a{z,e ie ) = y^a n (z)e ine with a n {z) = — [ a(z,6)e- in6 d6 (18) 
nez J 

Let f i — y f be the conjugation operator, determining the harmonic conjugate of a smooth func- 
tion. Defining the Bergman space if 2 as all complex-analytic and Lebesgue square-integrable 
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functions on the unit disc, then the orthogonal projection from L 2 (D + ) to H 2 is defined (e.g. 
|Gar0 9j) by the operator P via 

or explicitly 

nSZ n£Z+ 

The operator P : L 2 — > H 2 then can easily been seen as removing negative frequencies from 
the initial signal. 



7.2 Scaling Redux 

As usual we let \i[z, z) be real-analytic, absolute-summable and nonvanishing. Then 9* J2 P + q = n a vi zV ^ q 
takes the form 

el e fi = J2 c n(z)e in9 , with c r (z,z) = ^ a pq z p z« (19) 

q—p=r 

We conveniently now view the c n 's as Fourier coefficients of the function fi(z,6) i.e. Cj(z) = 
fi(z,j). Define the operator Pj-j on smooth functions via 



i 

= Y^c n (z)e ine 



n=k 



Let 1Z(Q) be the space of real-analytic functions of two variables on a region Q. Then define 
the following space; 



} \/z G fi/{0} there are infinitely many n such that 

C(Q) = <; g(z) G Tim ^ n)rg{z ^ + l) and limsup^ |fe±±)| < i 

Clearly 6{Q) is "most"© of K(Q) since it accounts for, in some sense, those real-analytic func- 
tions with "non-sparse" spectrums. The classical ratio test for infinite series ensures that \i G C 
are also absolute-summable. We have the following result about convergence on compact sub- 
sets. 

Proposition 7.1 Let K C VL be compact and let Q = [J, k « gZ2 Qk,i- Then Q{K) is dense in 
C(K) with respect to the uniform norm. 

Proof Let u)(z) G C{K) and let ^ Zj G supp u C K. Then, efu(z) = J2z u n{z)e me with 
oo n (zj) 7^ for infinitely many n G Z. We may pick an l(zj) such that u(zj,l(zj)) ^ 0. By 
the assumptions of C , there exists a finite k(zj) < —1 such that k(zj) + l(zj) + 2 > and 
< \u)(zj,k(zj))\ < \u(zj,l(zj))\, namely k(Zj) = —(l(zj) + 1). 



3 In the reasonable, informal way rather than a measure-theoretic sense 
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The varieties {z;u r (z,z) = 0} define (possibly degenerate) circles. Therefore, there is an 
e-neighborhood N e .(zj) = Nj around Zj on which there are two simple functions, — oo < kj(z) < 
k(zj) and l(zj) < lj(z) < oo on z E Nj, for which 

< < MzM^l ^ z e Nj/{0} 

and 

lj(z) + kj(z) + 2 > VzeNj 

Define 

Qj(z) = x(Nj)P kW{j) e*u 

where k(j) = min 2e jv J kj(z) and = max zg 7v i lj(z)- Then K d\jNj provides an open cover 
reducible to 



v 

supp u> C JVj 

3=1 

Consider the following function 

v 

3=1 

where ki(z) = mnij tZ kj(z) and lj(z) = max.j !Z lj(z) on z e C\Nj ^ in the case of overlapping 
neighborhoods. By design we have that 

with k{u) = min j kj(z) and 1{oj) = max j lj(z). Further, if Sr(z) is the R'th partial Fourier 
sum ofuu(z) notice that 

\Q p (z)-S R (z)\=0(±i), 5>0 
with L = R — max{\k(uj)\,l(uj)}. The classical bound 

\u(z) - S R (z)\ < ^ v n (z) 

\n\>R 

guarantees that on letting min{|/c|, 1} / oo and £j \ 0, that lim^oo Ct p (z) = ou(z) uniformly 
since the Fourier series can be brought as close as wanted in the mesh limit. 

The following corollary is then immediate. 

Corollary 7.2 Let K C f2 be compact and let H = U(fc i)ez 2 ^k,i- Then T-L(K) is dense in 
C(K) with respect to the uniform norm. 

7.3 Putting it all together 

Define 

0(Q) = {/iG 1Z(VL) satisfying condition 4 of H-ness} 

and let V = C(~)0. Then, if K is compact, for /x e T)(K) we see that X = VL p (z)d + Cl p (z)d can 
be chosen to approximate X = fxd+jid so that their integral curves are arbitrarily close in L P (K) 
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for 1 < p < oo via Poincare's inequality. Then | |x — x| |lp < C\\p — Q p (z)\\l p + | |(x — x)o| \u> < 
C'(e + diam{i^/supp/}). We then make the obvious choice setting K = supp/. 

Let p a fi be the seminorm on the Frechet space S of Schwarz-class functions on C, namely 
Pa,/?(0) — su Pxgc \x a d^(j)\, which generates the usual topology on S. If / G C%°(D + ) then clearly 
Igf(s) G S. We let s be the transverse flow induced by X ± from which the corresponding 
Hilbert transform H§ is defined. By continuity of Hilbert transforms on Schwartz functions, we 
assume that Uglf — H s lf f is small in the induced norm p on S. That being the case, then 
since Xg- = 9*a(z)*-^= for some function a(z) and because differentiation is continuous on S, 
we see that p{X e L (HsI e x — H s lf /)} may therefore be made as small as desired. 

We have therefore established the following corollary. 

Corollary 7.3 Let e > 0, / G C^°(D + ), p G V(D + ) and let Igf be the ray transform of f over 
the integral curves ofYg = 9*(pd + pd). Then there exists functions w e {z), and Cl p (z) such that 

\\I jc f-I X f\\L H D+)<Ce l<g<oo 

for C = C(supp f) constant and where Xg = 9 Jf — (Vl p d + Cl p d) . Suppose further that WHglg- f — 
HIef\\s < <K e )- Then there exists X\(z) such that 



sup 

z€D+ 



f(z) - ^ jf" ' P{\l,6)XiH{I e f){s{ze- ie U e )d6 



< e 



We may now summarize our stability and approximation results in the following theorem. 



Theorem 7.4 Lete > and H(Igf)(s(ze l0 ),e ie ) be given. Suppose thatp(z,z) = 
is a real- analytic function on D + and that f G C^°(D + ). Define c r (z,z) = Y2q- 



p+q>0 "'PI 



d"T)ri Z Z 



p?q 



q—p=r PI 



OrtnZ Z 



p?q 



Furthermore, suppose X*s, X*s z , and X^s z are meromorphic for X G D + . If, for all z G 
D + /{0}, there exist infinitely many j G Z such that Cj(z) and Cj +1 (z) are both nonzero, and if 



lim sup 



I C J + lO) I 

j^oo I Cj .( 2 ) 



< 1, then there exist functions w € (z), and Cl p (z) such that 
\\I*f -I X f\\mD+)<Ce l<g<oo 



If, in addition, we have \\HsI^ f — HIgf\\s < 5(e) then there exists a function Xi : z — > Xi(z) 
satisfying the following inequality 



< sup 

z&D+ 



/(*) 



w e (z) 
An 



2?r 



P(\i, 9)X e L H(I e f)(s(ze- w ), e w )d9 



< e 



(20) 



where Igf is the trace of f over the integral curves of Xg = 9*(pd + pd) and where Xq = 

i9±(-n p d + & p d). 
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